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ABSTRACT 

We report the results of a systematic timing analysis of all archival Rossi X-Ray 
Timing Explorer (RXTE) observations of the bright black-hole binary GRS 1915+105 
in order to detect high-frequency quasi-periodic oscillations (HFQPO). We produced 
power-density spectra in two energy bands and limited the analysis to the frequency 
range 30-1000 Hz. We found 51 peaks with a single trial significance larger than 3<r. 
As all but three have centroid frequencies that are distributed between 63 and 71 Hz, 
we consider most of them significant regardless of the number of trials involved. The 
average centroid frequency and FWHM are 67.3±2.0 Hz and 4.4±2.4 Hz respectively. 
Their fractional rms varies between 0.4% and 2% (total band detections) and between 
0.5% and 3% (hard ban detections). As GRS 1915+105 shows large variability on 
time scales longer than Is, we analysed the data in 16s intervals and found that the 
detections are limited to a specific region in the colour-colour diagram, corresponding 
to state B of the source, when the energy spectrum is dominated by a bright accretion 
disk component. However, the rms spectrum of the HFQPO is very hard and does not 
show a flattening up to 40 keV, where the fractional rms reaches 11%. We discuss our 
findings in terms of current proposed models and compare them with the results on 
other black-hole binaries and neutron-star binaries. 
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1 INTRODUCTION 

The large wealth of X-ray observations of black-hole tran- 
sients (BHTs) accumulated by the Rossi X-Ray Timing Ex- 
plorer (RXTE) mission has allowed us to obtain a much 
more precise view of the properties of accretion onto stellar- 
mass black holes (see e.g. Belloni 2010; Fender 2010; Belloni 
et al. 2012) compared to the sparse observations of a few ob- 
jects which were available before. In addition to detailed in- 
formation on low-frequency (0.1-30 Hz) Quasi- Periodic Os- 
cillations (QPO, see e.g. Motta et al. 2011), RXTE led to 
the discovery of higher-frequency features, which sample the 
range expected from signals associated to Keplerian motion 
in the innermost regions of an accretion disk around a black 
hole. These QPOs at higher frequencies (30-450 Hz) are 
known as High-Frequency QPOs (HFQPOs) and are mostly 
weak and elusive signals. Although the existing database of 
RXTE observations is very large, only a small number of de- 
tections is available, associated only to a few sources: GRS 
1915+105 (Morgan et al. 1997; Strohmayer 2001a; Belloni 
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et al. 2001; Remillard et al. 2002b; Belloni et al. 2006), GRO 
J1655-40 (Remillard et al. 1999; Strohmayer 2001b), XTE 
J1550-564 (Homan et al. 2001; Miller et al. 2001; Remillard 
et al. 2002a), H1743-322 (Homan et al. 2005; Remillard et 
al. 2006), XTE J1650-500 (Homan et al. 2003), 4U 1630-47 
(Klein- Wolt, Homan & van der Klis 2004), XTE J1859+226 
(Cui et al. 2000) and IGR J17091-3624 (Altamirano & Bel- 
loni 2012). Some of these detections appear in pairs. Re- 
cently, Belloni, Sanna & Mendez (2012) examined systemat- 
ically the full archive (with the exception of two sources, see 
below) and found only eleven statistically-significant detec- 
tions from two sources, after taking into account the number 
of trials. All detections appear when the sources are in an 
intermediate state and none of them together with a type-C 
LFQPO (see below). For all detections, there is an indication 
that there are preferred frequencies, but the small number 
of detections prevents a more solid assessment. 

Low-Frequency QPOs (LFQPOs) are much more com- 
mon in BHTs (see Motta et al. 2011 and references therein). 
They have been classified into three separate types (called 
A, B and C). Fype-C QPOs are the most commonly ob- 
served. They vary in frequency over a rather large range 
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and have been interpreted as related to the frame-dragging 
precession time scale (see e.g. Stella, Vietri & Morsink 1999; 
Ingram, Done & Fragile 2009). Type-B QPOs are restricted 
over a much small range of frequencies, while type-A QPOs 
are much fainter signals and less information is available. 

Two sources were not considered by Belloni, Sanna & 
Mendez (2012). ICR J17091-3624 and GRS 1915+105, as 
their peculiar behaviour requires a special analysis. IGR 
J17091-3624 was analysed in detail by Altamirano & Belloni 
(2012). Here we concentrate on GRS 1915+105 and present 
a systematic high-frequency timing analysis of the full set of 
existing RXTE observations of BHTs. 



2 GRS 1915+105 

GRS 1915+105, discovered in 1992 as a bright transient 
which has remained active since, is a bright black hole bi- 
nary which displays a very peculiar behaviour in the form 
of complex structured variability (see Belloni et al. 2000; 
Fender & Belloni 2004 for a review). Its luminosity is esti- 
mated to be near Eddington and it is the most important 
system to study the accretion-ejection connection. Its vari- 
ability was studied by Belloni et al. (2000), who classified 
its behaviour into 12 variability classes, extended later to 
14 (Klein- Wolt et al. 2002; Hannikainen et al. 2005). It was 
a unique source until in 2011 a new much fainter system 
was discovered, IGR J17091-3624, which displayed the very 
same type of peculiar variability (see Altamirano et al. 2012 
and references therein). 

As mentioned above, high-frequency oscillations were 
seen with RXTE from GRS 1915+105. Morgan et al. (1997) 
reported the discovery of a HFQPO around 65-67 Hz, which 
varied little between the very few observations where it was 
seen. Its total rms amplitude was low, around 1%, with a 
hard spectrum (up to 6% in the 15-25 keV band). Its quality 
factor Q, defined as the ratio between centroid frequency and 
FWHM was high, around 20. Belloni et al. 2001 analysed in 
detail two of these observations and found that the spectrum 
of the 67 Hz QPO changed dramatically as a function of a 
much slower 0.067 Hz oscillation, and discovered a broader 
(Q~3) QPO at 27 Hz, also dependent on the slower oscil- 
lation. Around the same time, Strohmayer (2001a) discov- 
ered a 41 Hz QPO when analyzing new observations taken 
in 1997 where a 69 Hz QPO was detected. In this case, 
the 41 Hz and 69 Hz peaks were detected simultaneously. 
It was noticed that the ratio between these frequencies ap- 
proached 5:3 (Kluzniak & Abramowicz 2002). Noticeably, 
the sequence 27:41:69 approaches the 2:3:5 ratios. Remil- 
lard et al. (2002b) reported the discovery of two additional 
peaks at 164 and 328 Hz (consistent with being harmonics), 
from a selection made in hardness and intensity from a set 
of RXTE observations in 1997. Belloni et al. (2006) reanal- 
ysed the data and confirmed the ~170 Hz , although with a 
rather low Q~2. 



3 DATA AND ANALYSIS 

We selected all RXTE observations of GRS 1915+105 avail- 
able in the archive from the start to the end of the mission, 
concentrating on the data from the Proportional Counter 



Array (PCA) instrument. We analyzed a total of 1807 ob- 
servations, for a total of 5181 ks exposure. Given the size of 
the available database, we followed a semi-automatic proce- 
dure similar to that followed by Belloni, Sanna & Mendez 
(2012): 

(i) We produced Power Density Spectra (PDS) from two 
energy bands: channels 0-79 and 14-79 (total and hard band 
respectively, corresponding to energies 1.9-29.8 keV and 5.5- 
29.8 keV for the first observations where we had a detection 
(see below) and to energies 2.1-33.4 keV and 6.1-33.4 keV 
at the end of the mission, due to gain changes in the detec- 
tors. For some observations, we could not use these channels 
bands due to incompatible data modes. In these cases, we 
chose the closest approximation to the desired bands. In 
a few cases, the data modes were such that we used the 
full channel range and/or the total and hard band were the 
same. The choice of the hard band was made in order to in- 
clude a sufficient number of photons and to allow the exam- 
ination of more observations (see next bullet). We consider 
a separate observation (or dataset) the data corresponding 
to a RXTE observation ID. 

(ii) The large number of observations made it necessary 
to perform an automatic search in order to select candidates 
for HFQPOs. For each observation, the adopted procedure, 
applied both to the total and hard energy band, was the 
following: 

• The dataset was subdivided into intervals of 16s du- 
ration, from each interval we produced a PDS and av- 
eraged all PDS. The time resolution for all spectra was 
4096 points per second, corresponding to a Nyquist fre- 
quency of ~2 kHz. The PDS were normalised according to 
Leahy et al. (1983) and rebinned in such a way that each 
frequency bin was larger than the previous one by ~2%. 
Uncertainties in power were estimated following van der 
Klis (1988). 

• Since our interest was only in the high-frequency re- 
gion, we limited our analysis to the frequency range 30- 
1000 Hz. Extending the threshold to lower frequencies 
(where the 27 Hz features was detected) would have cre- 
ated problems due to residual source noise, which would 
complicate the automatic detection. 

• We fitted the PDS with a model consisting of the 
sum of a power law (to account for the Poissonian noise 
component) and a Lorentzian. The initial procedure was 
completely automatic. The slope of the power law was 
not fixed to in order to fit any possible remaining tail of 
low-frequency features over the roughly flat Poisson level. 
We first fixed the Lorentzian centroid frequency to val- 
ues between 30 and 1000 Hz in steps of 1 Hz, limiting the 
FWHM between 0.5 Hz and 1000 Hz. Then, we made a fit 
with a free centroid around the frequency corresponding 
to the minimum chi square from the previous procedure 
(regardless of the quality of the fit), obtaining the best fit 
for each observation. This procedure is aimed at detect- 
ing a single HFQPO in each PDS and is unable to detect 
multiple peaks. A successful fit was obtained for all obser- 
vations at this stage, although for many cases the QPO 
detection was not significant. 

• We then estimated the significance of the detection 
by dividing the normalisation of the Lorentzian (corre- 
sponding to the integral over positive frequencies) by its 
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la error on the negative side. We retained only the 109 
detections where this significance was larger than 3. 

• All the 109 detections were examined visually, since 
in many cases the fits did not make sense. For exam- 
ple, in some cases the FWHM resulted in a very broad 
or very narrow feature (hitting the imposed limits), or 
the centroid frequency itself was also close to the search 
limits, or finally there was remaining source noise above 
30 Hz, leading to a spurious QPO detection. In addition, 
we discarded all fits where the QPO had a quality factor 
Q — vq/FWHM significantly less than 2. 

• For the observations that survived the procedure, we 
repeated the fit manually and we estimated the signifi- 
cance of the fit. In a large number of cases, we found a 
peak at a consistent frequency both in the hard and total 
PDS. In these cases, we retained the parameters of the 
fit with the highest significance. No pairs of peaks at in- 
consistent frequencies between total and hard band were 
found. We were left with 51 observations containing a sig- 
nificant peak. For all our final fits, the reduced chi square 
was close to unity. The log of observations is shown in 
Tab. [T]and the HFQPO parameters in Tab. H 



(iii) The final HFQPOs have a significance larger than 3cr, 
but for a single trial. To keep into account the number of tri- 
als in our procedure, we had to consider the number of obser- 
vations analysed and the number of independent frequencies 
in the PDS. As number of observations, we used 1807 if the 
QPO was detected in both bands, as the second detection 
did not add significance to the first. If the QPO was detected 
in only one of the two bands, we used 1807*2. To estimate 
of the number of independent frequencies scanned for each 
observation where a QPO was detected, we divided the fre- 
quency range used for the analysis (970 Hz) by the FWHM 
of the QPO detected in that particular observation. Follow- 
ing this procedure, the number of trials results very large, of 
the order of 10 6 , which lowers considerably our sensitivity. 
Only 4 detections remain with a final detection probability 
less than 1%. However, 49 of our 51 detections fall in the fre- 
quency range 58.9-71.3 Hz and 48 of them in the frequency 
range 63.5-71.3 Hz, suggesting that most of them are real. 
Using that reduced frequency range (14 Hz) increases our 
<1% observations to 14. We report all observations with 
single trial > 3<r as a posteriori the chance of having al- 
most the totality of the QPOs within such a small range 
of frequencies is rather low. Two detections are at higher 
frequency >100 Hz and have a high probability of being 
spurious after number of trials are taken into account. We 
include them in the following discussion not having strong 
statistical reasons to discard only these two values (we note 
that their frequency is consistent with being twice some val- 
ues obtained in other observations, 67.05 Hz and 71.49 Hz). 
We checked whether these two observations also feature a 
lower-frequency HFQPO, but found none. 

(iv) Notice that our procedure can only detect one peak 
per PDS. In principle, two different features can be detected 
in the total and hard PDS, but this did not happen in any 
of our observations. 



Table 1. List of GRS 1915+105 observations with a HFQPO de- 
tection. Columns are: observation #, observation ID, observation 
date (MJD), variability class, exposure time 



Obs. N 


ObsID 


MJD 


Class 


Exp. (s 


1 


10408-01- 


01-01 


50179 


.236 


K 


5600 


2 


10408-01- 


-03-00 


50190 


.578 


7 


4768 


3 


10408-01- 


-04-00 


50193 


.439 


7 


8768 


4 


10408-01- 


-05-00 


50202 


.845 


7 


9584 


5 


10408-01- 


-06-00 


50208 


.584 


7 


9856 


6 


10408-01- 


-07-00 


50217 


.656 


7 


9984 


7 


10258-01- 


-10-00 


50351 


.668 


1-' 


5488 


8 


20186-03- 


-01-03 


50628 


.615 


K 


3504 


9 


20186-03- 


-01-04 


50629 


.483 


K 


10608 


10 


20402-01- 


-38-00 


50649 


.426 


7 


7472 


11 


20402-01- 


-39-00 


50654 


.028 


7 


6864 


12 


20402-01- 


-39-02 


50658 


.511 


7 


2128 


13 


20402-01- 


•41-02 


50679 


.372 


S 


3312 


14 


20402-01- 


-54-00 


50763 


.209 


8 


10352 


15 


20402-01- 


-55-00 


50769 


.229 


8 


8768 


16 


20402-01- 


-56-00 


50774 


.224 


8 


9408 


17 


20402-01- 


-57-01 


50786 


.362 




4672 


18 


20402-01- 


-58-01 


50788 


.366 


H 


4512 


19 


20402-01- 


-60-00 


50804 


.908 


8 


12368 


20 


30703-01- 


-01-00 


50810 


.027 


8 


5168 


21 


30402-01- 


-01-00 


50820 


.164 


8 


6880 


22 


40403-01- 


-06-00 


51158 


.274 


K 


2144 


23 


40703-01- 


-12-00 


51288 


.149 


U! 


9344 


24 


40403-01- 


-07-00 


51291 


.081 


W 


2816 


25 


40703-01- 


-13-01 


51299 


.134 


U! 


6080 


2(5 


40703-01- 


-26-00 


51407 


.769 


U! 


8320 


27 


40703-01- 


-28-00 


51418 


.781 


W 


1536 


28 


40703-01- 


-28-02 


51418 


.920 


W 


1856 


29 


40703-01- 


-29-00 


51426 


.769 


U! 


1872 


30 


40703-01- 


-29-01 


51426 


.839 


U! 


1504 


31 


40703-01- 


-29-02 


51426 


.909 


U! 


1920 


32 


40703-01- 


-31-00 


51447 


.728 


8 


9984 


33 


50703-01- 


-10-02 


51681 


.213 


P 


864 


34 


50703-01- 


-60-03 


52052 


.234 


P 


1040 


35 


50703-01- 


-61-01 


52059 


.119 


P 


1648 


36 


60701-01- 


-09-00 


52192 


.641 


V 


6976 


37 


70702-01- 


-55-00 


53361 


.525 


5 


4832 


38 


80701-01- 


-28-00 


52933 


.627 


6 


1952 


39 


80701-01- 


-28-01 


52933 


.696 


8 


1616 


40 


80701-01- 


-28-02 


52933 


.764 


8 


1008 


41 


80127-04- 


-02-00 


52945 


.238 


8 


15104 


42 


80701-01- 


-31-00 


52955 


.570 


6 


11440 


43 


90105-08- 


-02-00 


53503 


.870 


8 


5104 


41 


90701-01- 


-38-00 


53347 


.564 


8 


1952 


45 


91701-01- 


-12-00 


53515 


.075 


8 


2112 


46 


91701-01- 


-12-01 


53515 


.140 


8 


2352 


47 


92092-03- 


•01-00 


53707 


.683 


8 


13088 


48 


93701-01- 


-01-00 


54285 


.950 


7 


1648 


49 


93701-01- 


-06-00 


54320 


.822 


8 


3056 


50 


93411-01- 


-01-00 


54326 


.709 


8 


3424 


51 


95701-01- 


-20-00 


55335 


.140 


P 


1216 



3.1 Colour— colour and hardness— intensity 
diagrams 

We use the 16-s time-resolution Standard 2 mode data to 
calculate X-ray colours. For each of the five PCA detectors 
(PCUs) we extracted count rates every 16 seconds in the 2.0- 
6.0 keV, 6.0-16.0 keV, 16.0-20.0 keV and 2.0-20 keV bands. 
To obtain the count rates in the exact energy ranges, we in- 
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Table 2. List of detected HFQPOs for GRS 1915+105. Columns are: observation number, bands in which the QPO was 
detected (Total/Hard: in boldface the most significant, from which the parameters were obtained), centroid frequency, 
FWHM and fractional rms of the QPO, chance probabilities for single trial (with number of sigma), full frequency range 
and restricted frequency range (see text), net source average observed count rate (in PCU2, with average number of active 
PCUs. Probabilities of and 1 indicate values too close to those values to be displayed. 



Obs N 



B 



v (Hz) FWHM (Hz) 



%rms 



Po 



Rate (Npcu) 



1 


1 II 


66 


21 


+0.53 
-0.75 


3.11 


+1 
-2 


60 
11 





58 


+0 
-0 


10 
12 


2 


1 II 


68 


79 


+0.40 
-0.68 


2.86 


+ 2 
-2 


11 
50 





62 


+0 
-0 


0!) 
14 


3 


1 II 


67 


34 


+0.14 
-0.12 


3.47 


+ 
-0 


44 
50 





86 


+0 
-0 


04 
04 


4 


TH 


66 


31 


+0.31 
-0.29 


3.15 


+0 
-0 


54 
75 





95 


+0 
-0 


07 
(19 


5 


TH 


65 


10 


+0.13 
-0.13 


4.08 


+0 
-0 


29 
38 


1 


21 


+0 
-0 


03 
05 


6 


T 


66 


98 


+ 
-0 


fil 

50 


2.66 


+1 
-1 


04 

53 





56 


+0 
-0 


09 
09 


7 


TH 


67 


71 


+0 
-0 


58 
59 


5.16 


+1 

-2 


31 
07 





41 


+0 
-0 


07 
08 


8 


H 


63 


54 


+1 

-1 


38 
42 


9.48 


+3 
-6 


65 
28 


1 


17 


+0 
-0 


20 
29 


9 


TH 


68 


43 


+0.89 
-0.88 


10.63±3 


42 
43 


1 


06 


+0 
-0 


11 
14 


10 


TH 


67 


40 


+0.59 
-0.60 


3.52 


+ 1 

-2 


67 
20 





42 


+0 
-0 


00 
08 


11 


I'll 
l ±1 


68 


92 


+0.21 
-0.21 


4.20 


+ 
-0 


68 
80 


1 


13 


+0 
-0 


00 
07 


12 


TH 


68 


92 


+0.39 
-0.41 


3.60 


+1 
-1 


24 

43 





98 


+0 
-0 


11 
13 


13 


TH 


65 


79 


+0.41 
-0.49 


5.13 


+1 

-1 


23 

44 





96 


+0 
-0 


08 
09 


14 


TH 


65 


96 


+ 1.17 
-1.02 


12.03+3 


82 
71 





91 


+0 
-0 


09 
11 


15 


TH 


68 


07 


+0.25 
-0.24 


3.99 


+0 
-0 


59 
73 


1 


01 


+0 
-0 


00 
07 


16 


TH 


68 


72 


+0.37 
-0.38 


4.90 


+1 

-1 


26 
91 





82 


+0 
-0 


08 
10 


17 


H 


62 


65 


+ 1 
-1 


40 
35 


10.42+J 


39 
98 


1 


06 


+0 
-0 


12 
13 


18 


H 


58 


90 


+1 

-1 


67 
80 


17.62+j 


91 

87 


1 


35 


+0 
-0 


15 
16 


19 


TH 


64 


75 


+0.31 
-0.33 
+0.51 
-0.53 


3 


91 


+ 

-1 


86 
11 





62 


+0 
-0 


05 
ill. 


20 


TH 


65 


89 


3 


87 


+1 

-2 


48 
14 





50 


+0 
-0 


07 
09 
05 
05 


21 


TH 


64 


99 


+0.24 
-0.24 


3 


21 


+ 
-0 


64 
76 





66 


+0 
-0 


22 


T 


65 


45 


+0.82 
-0.95 


2 


55 


+1 
-1 


75 
95 





79 


+0 
-0 


13 

13 


23 


TH 


70 


19 


+0.49 
-0.42 


3 


76 


+1 
-1 


36 
80 


1 


03 


+0 
-0 


13 
15 


24 


H 


70 


39 


+0.49 
-0.48 


2 


97 


+ 2 
-1 


10 
94 


1 


17 


+0 
-0 


18 
22 


25 


H 


68 


92 


+0.57 
-0.41 


3 


66 


+ 1 
-2 


49 
05 





91 


+0 
-0 


13 
15 


26 


TH 


68 


89 


+0.56 
-0.67 


5 


41 


+ 1 
-1 


31 
08 


1 


39 


+0 
-0 


14 
17 


27 


TH 


68 


82 


+0.39 
-0.39 


3 


82 


+1 

-1 


17 

38 


1 


66 


+0 
-0 


17 
19 


28 


H 


70 


26 


+0.45 
-0.45 


2 


98 


+1 

-1 


00 
41 


1 


32 


+0 
-0 


18 
21 


29 


H 


68 


57 


+ 1.40 
-1.40 


7 


04 


+ 2 
-4 


87 
34 


1 


40 


+0 
-0 


22 
22 


30 


H 


69 


58 


+0.47 
-0.54 


3 


01 


+1 
-1 


42 
97 


1 


29 


+0 
-0 


19 
26 


31 


TH 


69 


31 


+0.34 
-0.32 


2 


95 


+ 
-1 


86 
09 


1 


44 


+0 
-0 


15 
16 


32 


TH 


69 


61 


+ 1.08 
-0.76 


6 


60 


+ 2 
-4 


31 

34 


1 


04 


+0 
-0 


15 
21 


33 


T 


134.10 


+ 1 
-1 


.88 
.80 


6 


81 


+ 2 
-2 


20 

96 


2 


00 


+0 
-0 


30 
33 


34 


TH 


67.99 


+1 

-1 


10 
04 


5 


28 


+1 
-3 


97 
78 


1 


63 


+0 
-0 


21 
21 




TH 


142.98 


+3.81 
-3.15 


18.71 
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Figure 1. Top: Colour-colour diagram for all RXTE/PCA ob- 
servations of GRS 1915+105. Each point is accumulated for 16 
seconds. The magenta and black circles mark all segments in- 
cluded in the observations where a HFQPO was detected. Dif- 
ferent colours mark point in different areas (see text). Bottom: 
corresponding hardness-intensity diagram. 



terpolated linearly between PCU channels. Each light curve 
was corrected for dead-time effects following the methods 
suggested by the RXTE team. We also remove instrumen- 
tal drop-outs and subtracted the background contribution 
in each band using the standard bright source background 
model for the PCA, version 3.c£J We defined the soft colour 
(HR1) as the count rate in the 6.0-16.0 keV band divided by 
the rate in the 2.0-6.0 keV band and the hard colour (HR2) 
as the ratio of the count rates in the 16.0-20.0 keV rate di- 
vided by the 2.0-6.0 keV rate. We used the count rates in the 
2.0-20 keV band as a measurement of the source intensity. 

In the 16 years of activity, the RXTE gain changed 
5 times with each new high voltage setting of the PCUs 
(Jahoda et al. 2006). In order to correct for these gain 
changes as well as the differences in effective area between 



1 PCA Digest at http://heasarc.gsfc.nasa.gov/ for details of the 
model 



the PCUs, we used the method introduced by Kuulkers et 
al. (1994): for each PCU we calculate, in the same manner 
as for GRS 1915+105, the colours of the Crab, which can 
be supposed to be constant. We then average the 16s Crab 
colours and intensity for each PCU for each day. For each 
PCU we divide the 16s colour and intensity values obtained 
for GRS 1915+105 by the corresponding average Crab values 
that are closest in time but in the same RXTE gain epoch. 
Then, we average the colours and intensity over all active 
PCUs to obtain colours and intensities every 16 seconds. 

The colour-colour diagram (CCD) and hardness- 
intensity diagram (HID) obtained with the 16s points are 
shown in Fig. [1] The colour dots marks the 16s segments be- 
longing to observations where we detected a HFQPO. The 
colours are assigned in the following way: 

• magenta: (HRi > 0.7 & HR 2 < 0.32) or (HRi > 0.9 & 
HR 2 < 0.4) 

• black: remaining points. 



4 RESULTS 

4.1 Where HFQPOs are detected 

Figure [2] shows the distribution of the parameters of our 
detected peaks. The centroid frequencies are concentrated 
around 67 Hz, with 48 our of 51 detections between 63.5 
and 71.3 Hz, with an average of 67.3+2.0 Hz. The FWHM 
peaks between 3 and 4 Hz, with an average of 4.4+2.4 Hz. 
The fractional rms, which is of course affected by our sensi- 
tivity in its lowest values, peaks at 1 %. In Fig. [3] we show 
two examples of detections: that with the highest single-trial 
significance (Obs. 34, no- = 18.5) and one of the 3a ones (Obs. 
1, n CT =3). 

Our 51 detections listed in Tab. [2] come from observa- 
tions from seven variability classes: k,^,ijl,5,oj,p,v (see Belloni 
et al. 2000; Klein- Wolt et al. 2002). For many of these, the 
average rates quoted in Tab. [2] are only indicative, as there 
are large rate variations during the observation. Noticeably, 
none of the detections corresponds to observations in class 
X, which is the one showing consistently a type-C QPO. This 
is in agreement with the results of Belloni, Sanna & Mendez 
(2012), who found no HFQPO and type-C detection in other 
sources. Belloni et al. (2000) have shown that the >ls vari- 
ability of GRS 1915+105 can be reduced to the transition 
between three basic spectral states, called A, B and C. Only 
two variability classes do not involve state transitions: <f> and 
X, none of which appear in our sample. 

In order to establish whether the presence of the 
HFQPO is related to a specific state, we extracted PDS in 
the total band from selections made on the HID and CCD 
in Fig. [1] All 16s segments from the observations where we 
detect a HFQPO are shown in colour in Fig.[T] We identified 
two regions (for their precise definition see above): the soft 
points (magenta), which occupy a dense region of the CCD, 
and the remaining points (black) . 

We extracted the average PDS from the two regions 
from both the total and hard bands, shown in Fig. [4] It 
is clear that the HFQPO is detected only in the magenta 
region, which corresponds to state B (the CCD in Fig. [T] 
is comparable with that in Belloni et al. 2000, with axes 
flipped). The black points cover both a softer region of state 



6 T. M. Belloni & D. Altamirano 



10 



s 

"S 

2 4 



50 60 70 



80 90 100 110 120 130 140 150 
HFQPO Frequency (Hz) 



10 



2 5 



nnn n 



nn 



6 8 10 12 14 
HFQPO FWHM (Hz) 



16 18 20 



~ 4 



o 



0.5 



1.0 



1.5 2.0 
HFQPO % rms 



2.5 



3.0 



3.5 



Figure 2. Distributions of best-fit parameters for all HFQPOs 
detected (see Tab. From top to bottom: centroid frequency, 
FWHM and fractional percentage rms. 



B and part of state C. No state A point is either black or 
magenta. The magenta selection yields an averaged QPO 
centered at 66.9T0T Hz, with a FWHM of 5.5T0.3 Hz and 
a fractional rms of 0.69±0.02% for the full band (signifi- 
cance 24.4<t) and centered at 67.1T0T Hz, with a FWHM 
of 5.9T0.3 Hz and a fractional rms of 0.93±0.02% for the 
hard band (significance 26.3<r). The black selection, assum- 
ing the same shape for the QPO, yields 3a upper limits of 
0.8% (total band) and 1.0% (hard band). 

Having identified the region in the CCD where HFQ- 
POs are found, we accumulated a PDS from all 16s points 
in that region from observations where we did not detect a 
HFQPO in single observations. This selection resulted in a 
positive 7.9ct detection (see Fig. [5}, indicating that state B is 
indeed associated to the presence of HFQPO. The centroid 
is 65.8±0.3 Hz, the FWHM 4.3±0.6 Hz and the integrated 
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Figure 3. Plot of our most significant (top panel, Obs. #39) 
and one of the least significant (bottom panel, Obs. #1) HFQPO 
detections. 



fractional rms 0.33±0.02%, lower by roughly a factor of three 
than that from the magenta points. It corresponds to a sig- 
nificance of 8.3ct. Not all PDS averaged need to contain an 
identical signal and therefore this value had to be taken as 
an average value. If there are PDS with a weaker signal, 
this value must be seen as an upper limit. Notice that also 
the centroid frequency is significantly lower, although it is 
difficult to say anything more as these are averages over a 
large number of observations. The same result was obtained 
by Altamirano & Belloni (2011) for ICR J10791-3624. The 
HQFPO in that source was found only in a subset of bright 
observations; selecting all remaining bright observations did 
reveal a similar QPO peak. 

In order to check for a possible correlation between 
QPO frequency and rate, we averaged the Crab-corrected 
rates from the magenta points observation by observation, 
obtaining a clean count rate for each observation, limited to 
the intervals where the HFQPO is detected. A linear regres- 
sion finds a correlation factor of -0.29, indicating the absence 
of a statistically significant correlation. 



4.2 Rms spectrum of the HFQPO 

In order to extract the rms spectrum of the QPO we selected 
two groups of observations. The first is single observation 
#5 (from 1996), the same analyzed by Morgan et al. (1997). 
The second is the combination of observations #38, #39, #40 
(from 2003), which are from the same day and all show a 
very strong HFQPO peak at a very similar frequency. The 
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Figure 4. PDS averaged over the colour regions in Fig. [JJ (corresponding to the observations where a HFQPO is detected) for the total 
and hard band (left and right, respectively). The insets in the top panels shows a zoom around the QPO area. No QPOs are detected in 
the bottom panels. Error bars in the top panels are smaller than the symbols. 



1.98 



"Non-detections" 



100 

Frequency (Hz) 



Figure 5. PDS obtained from the average of all 16s points in the 
magenta region, but belonging to observations where we do not 
detect a HFQPO in the single full observations. 
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Figure 6. Rms spectra of the HFQPO for observation #5 (black 
dots, 1996, see also Morgan et al. 1997) and for the combined 
observations #38,#39,#40 (gray stars, 2003). 



resulting rms spectra are shown in Fig. [6] The rms vs. energy 
dependence is very similar and with the 2003 observations, 
but we are able to extend it to 40 keV, showing that it 



continues to increase up to at least those energies, with an 
rms larger than 10 %. 
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4.3 The 41 Hz QPO 

Strohmayer (2001a) reported the discovery of a second QPO 
peak around 41 Hz for a set of five observations in 1997. We 
do not see this peak in single observations, although in each 
of them we detect a ~67 Hz HFQPO. However, averaging 
those five observations we reproduce the 41 Hz detection. 
These observations cover a subset of colours of the magenta 
area. Selecting all other points in this subset does not reveal 
the oscillation. The conclusion is that the 41 Hz feature is 
a much more transient phenomenon, unlike the 67 Hz one, 
which appears much more frequently. 



5 DISCUSSION 

We have performed a full analysis of 1807 RXTE/PCA ob- 
servations of GRS 1915+105 with the aim of obtaining a 
homogeneous set of detections of its HFQPOs. Our proce- 
dure returned 51 observations with a single-trial detection 
significance larger than 3a. All but three of these detec- 
tions correspond to centroid frequencies between 63.5 and 
71.3 Hz, one is at a lower frequency (59 Hz) and two are at 
higher frequencies (127 and 143 Hz). The latter are broadly 
consistent with being a factor of two higher in frequency 
than the ~67 Hz main group. The 48 detections are grouped 
around a centroid frequency of ~67 Hz, a FWHM of ^3 Hz 
and have a total fractional rms of 0.5-1.5%. 67 Hz appears 
to be a very special frequency associated to this particu- 
lar source. An analysis made selecting specific intensities 
and X-ray colours indicate that the HFQPO is associated 
only to state B of GRS 1915+105 (see Belloni et al. 2000). 
Three separate states have been defined for the variability 
of GRS 1915+105 (Belloni et al. 2000). State C is equivalent 
to the Hard Intermediate State of more conventional black- 
hole transients, where the spectrum is composed by a ther- 
mal disk component plus a relatively strong hard tail, which 
can dominate the flux (see Belloni 2010; Belloni, Motta & 
Munoz-Darias 2012). States A and B correspond to an en- 
ergy spectrum strongly dominated by the thermal disk com- 
ponent and have a PDS almost featureless (Reig, Belloni & 
van der Klis 2003) and can be compared to the "anomalous" 
state of black-hole binaries (Belloni 2010). Here the energy 
spectrum is dominated by a bright thermal component and 
the inner radius of the accretion disk is close to its lower 
value (Belloni et al. 1997a,b; 2000). We obtain no detection 
during state C, when strong type-C QPOs are observed (see 
Markwardt, Swank & Taam 1999; Fender & Belloni 2004 
and references therein). However, the rms spectrum of the 
HFQPO is very hard and increases up to 40 keV without 
any sign of flattening (see also Morgan et al. 1997). 

It is remarkable that the frequency we find is very 
close to that discovered from the only known X-ray source 
which displays properties similar to those of GRS 1915+105, 
IGR J17091-3624 (Altamirano & Belloni 2012). The fre- 
quency recovered when analysing all magenta points in 
GRS 1915+105 is 67.1+0.1 Hz, while the detection in IGR 
J17091-3624 is at 66.5+0.5 Hz, compatible within la. Al- 
tamirano et al. (2012) suggested that IGR J17091-3624 has 
a lower mass than GRS 1915+105, but even assuming a 
similar mass the detection of precisely the same frequency 
is puzzling. 



RXTE observed dozens of black-hole transients in its 
sixteen years of operation. Overall, the number of detections 
of HFQPOs from these systems was very limited (Belloni, 
Sanna & Mendez 2012). GRS 1915+105 alone yielded more 
than 40 detections, which seems to deviate from the pat- 
tern. However, it was shown that the detectability of HFQ- 
POs strongly depends on source state. All detections in other 
transients were associated to relatively rare states at high lu- 
minosity (the Soft-Intermediate and Anomalous states, see 
Belloni 2010). GRS 1915+105 does not follow the same pat- 
tern as conventional transients, but it is extremely bright 
and it was shown that its properties can be associated to 
the same basic states (see e.g. Reig et al. 2003; Soleri et al. 
2008; Belloni 2010). Although it is difficult to make a quan- 
titative comparison, it is clear that GRS 1915+105 spends 
a much higher percentage of time in HFQPO-related states 
and therefore its properties might be compatible with those 
of the others. 

The HFQPO frequency cannot be the Keplerian fre- 
quency at the Innermost Stable Orbit (ISCO) around the 
black hole in GRS 1915+105. For a dynamical mass of 
14.0+4.4 M (Harlaftis & Greiner 2004), the lowest pos- 
sible Keplerian frequency at the ISCO, obtained in the case 
of M=18.8 Mq and no spin, is 110 Hz. In order to reach 67 
Hz, the black hole mass should be as high as 30 Mq. How- 
ever, theoretical models proposed for the HFQPOs do not 
always associate the frequency to the ISCO. The relativistic 
resonance model (Kluzniak & Abramowicz 2002) associates 
the QPO frequencies to relativistic time scales at a specific 
radius, where these frequencies are in resonance, resulting in 
special frequency ratios. We found only single peaks in our 
analysis, but the fact that the variations in this frequency 
over 16 years of RXTE life was very small is qualitatively 
compatible with that idea. From Fig. [2] the QPO frequency 
is distributed over a narrow range of frequencies, but even 
limiting ourselves to the main peak in the distribution it 
varies ±4 Hz, i.e. 6% These variations must be evaluated 
in the framework of the model. The relativistic precession 
model (see Stella, Vietri & Morsink 1999 and references 
therein) makes predictions about the frequencies of both 
high-frequency peaks. Having only a single peak does not 
allow any testing of this model. The same applies to other 
models (disko-seismic: Nowak & Wagoner 1991; Accretion- 
Ejection Instability: Tagger & Varniere 2006 ; Inner- Torus 
Oscillation: Rezzolla et al. 2003). 

There is a detection of a high-frequency feature at a 
higher frequency in GRS 1915+105, 170 Hz (Remillard et 
al. 2002b, Belloni et al. 2006), although with a low qual- 
ity factor. This feature was found by adding observations 
belonging to variability class 9 and selecting only the hard 
intervals, corresponding to state C. This is very different 
from all other detections, in fact class 8 does not even fea- 
ture state-B intervals, but it is made by oscillations between 
states A and C (Belloni et al. 2000). This also means that the 
170 Hz feature is simultaneous with a type-C QPO, always 
present in state C. 170 Hz is high enough to be a Keple- 
rian frequency at the ISCO around the black hole in GRS 
1915+105, but it is clearly not the main feature observed in 
this source. 

Belloni, Sanna & Mendez (2012) compared the distribu- 
tion of HFQPO detected in black-hole transients with that 
of kHz QPOs in neutron-star binaries, using as compari- 
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son the large number of detections obtained for 4U 1636-63 
(Sanna et al. 2012). From the small number of detections, 
we cannot exclude that the two samples come from a sim- 
ilar distribution. In our case, the number of detections is 
larger and concentrated in a narrow range around 67 Hz. 
Comparing the histogram in Fig. [2] with their Fig. 5, the 
two distributions appear to be markedly different, although 
a number of selection effects could be at work. 

In 2013, the indian X-ray astronomy satellite AS- 
TROSAT will be in operation and will pick up RXTE's 
legacy on fast timing analysis (Agrawal 2006). In partic- 
ular, the higher effective area above 20 keV of the Large 
Area Xenon Proportional Counter (LAXPC) will be ideal for 
HFQPOs (see Fig.©. GRS 1915+105 in unpredictable in its 
behaviour, but hardness/intensity selection according to our 
data will allow to limit the analysis to the most promising 
intervals. 

In the more distant future, if selected by ESA the LOFT 
satellite will allow us to sample much fainter signals and will 
most likely lead to many more detections. Depending on the 
model, different frequencies are expected. However, it is clear 
that the strongest signal in GRS 1915+105 is around 67 Hz, 
which appears to be a special value for this source. Overall, 
the results of the complete analysis of HFQPOs in the BHTs 
observed by RXTE show that these signals are detectd only 
during specific source states. An observation strategy aimed 
at these states will allow to maximise detections. 
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